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Abstract. We study correlation functions of various static meson operators of size r at non-zero tempera- 
ture in 2+1 flavor QCD, including Coulomb gauge fixed operators and Wilson loops with smeared spatial 
parts. The numerical calculations are performed on 24 3 x 6 lattices using highly improved staggered quark 
action. We discuss possible implications of our findings on the temperature dependence of the static energy 
of QQ pair. 



1 Introduction 

At high temperature strongly interacting matter under- 
goes a transition to a new state called quark gluon plasma 
(QGP). Creating and exploring such state of matter is a 
subject of a large experimental program, see e.g. Ref. pQ. 
At the same time properties of strongly interacting matter 
at high temperatures can be studied in ab-initio calcula- 
tions that rely on lattice regularization of QCD (see Ref. 
[2J[3] for recent reviews). 

The suppression of quarkonium production has been 
suggested by Matsui and Satz as probe of QGP forma- 
tion in heavy ion collisions [4 . Interaction between heavy 
quark and anti-quark is expected to be modified by the 
deconfined medium at high temperature eventually lead- 
ing to the dissociation of the quarkonium states. The sup- 
pression of the quarkonium yields in heavy ion collisions 
was indeed observed experimentally (see Ref. [5]). The 
interpretation of the experimental findings, however, re- 
quires the knowledge of quarkonium properties at high 
temperatures among other things (see Refs. [SHU] for re- 
cent reviews). Therefore there has been a significant ef- 
fort to study in-medium properties of heavy quarkonium 
in recent years including lattice QCD studies. In-medium 
quarkonium properties are encoded in the spectral func- 
tions. A commonly used approach to obtain quarkonium 
spectral functions relies on lattice QCD calculation of the 
correlation functions in Euclidean time direction and ex- 
traction of the quarkonium spectral functions using the 
Maximum Entropy Method pTUlinill2llT3TTl] . However, 
it turns out that temporal quarkonium correlators are re- 
markably insensitive to the in-medium modification of the 
corresponding spectral functions due to the limited ex- 
tent of the imaginary time direction r < 1/ (2T) |T5] . The 
dissolution of the bound state peaks in the spectral func- 
tions is compensated by large enhancement of the spectral 
function in the threshold region in such a way that cor- 



responding Euclidean correlation functions do not change 
significantly for r < 1/ (2T) [TB] . This picture is corrobo- 
rated by the study of the spatial charmonium correlation 
functions 17 . Unlike temporal correlators spatial corre- 
lators can be calculated for separation larger than 1 / (2T) 
and show strong temperature dependence consistent with 
dissolution of the bound states [17] . 

Because of the above problems the calculation of the 
quarkonium spectral functions from finite temperature Eu- 
clidean time correlation functions is very difficult. An al- 
ternative approach is to calculate quarkonium spectral 
functions relies on effective theories. Quarkonium is char- 
acterized by three different energy scales: the heavy quark 
mass m, the inverse quarkonium size 1/r ~ mv and the 
quarkonium binding energy mv 2 . Integrating out the largest 
energy scale ~ m leads to an effective theory called non- 
relativistic QCD (NRQCD). This effective theory was com- 
bined with lattice QCD method to study quarkonium prop- 
erties both at zero [T8] and non-zero temperatures [TH 
[20] . Integrating out the scale mv gives an effective theory 
called potential non-relativistic QCD (pNRQCD) [211122] . 
The Lagrangian of pNRQCD is formulated in terms of sin- 
glet and octet fields of the heavy quark anti-quark pair and 
the heavy quark potential enters as the parameter of the 
effective Lagrangian. At non-zero temperature additional 
scales T and gT appear. Depending on how these scales 
are related to the above energy scales the heavy quark 
potential may be modified by the medium accordingly 
[23] . At non-zero temperature the potentials also have an 
imaginary part [231124] . The above effective field theory 
approach is based on the weak coupling, namely on the 
assumption that mv, T, gT 3> Aqcd- Within this frame- 
work in-medium quarkonium properties can be calculated 
1251 and one can also take the static limit and calculate 
the energy of a static QQ pair which like the potential 
also has an imaginary part. 
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It is not clear if in the interesting temperature range 
the above scale separation holds. However, we may ex- 
pect that the binding energy of quarkonium states is re- 
duced with increasing temperature. Close to the disso- 
lution point the binding energy is the smallest scale in 
the problem and all other scales, including Aqcd can be 
integrated out. In that regime the potential of the corre- 
sponding effective theory is equal to the energy of QQ pair 
[2"5] . Therefore calculation of the static energy at non-zero 
temperature in the non-perturbative region could be very 
helpful for the determination of in-medium quarkonium 
properties. 

At zero temperature the energy of static QQ pair is 
calculated by studying correlation function of static me- 
son operators at large Euclidean time separations. This is 
not possible and non-zero temperature due to the limited 
temporal extent r < 1/T. It has been suggested to extract 
the static quark potential at non-zero temperature from 
the spectral representation of the simplest finite tempera- 
ture static meson correlators, namely rectangular Wilson 
loops and MEM [271128]. The analysis was performed in 
quenched approximation. In this paper we are going to 
analyze the temperature dependence of correlation func- 
tions of various static meson operators in 2+1 flavor QCD 
across the chiral crossover region and discuss its implica- 
tion for the temperature dependence of the static QQ en- 
ergy. Apart from quantitative determination of the static 
energy at non-zero temperature the study of static meson 
correlators is interesting because they are simpler than 
quarkonium correlators and can provide some insight into 
the inter-quark interaction in the medium and the weakly 
interacting nature of QGP. One of the aims of the present 
study is to find out to what extent the temperature de- 
pendence is influenced by the choice of the interpolating 
meson operator. Obviously only the features that are not 
sensitive to the specific choice of the meson operators are 
of some relevance. Some preliminary results have been re- 
ported in conference proceedings [29,30 . 



2 Static quark anti-quark correlators in lattice 
QCD 

Let us consider a static meson operator of the form 

0(x,y;r) = ip\x,T)U(x,y,T)tp(y,T). (1) 

Here ip and tfr are the static quark and anti-quark fields, 
and U(x, y; r) is the spatial Wilson line connecting points 
a; and y. After integrating out the static quark anti-quark 
fields the correlation function of this operator at time r 
becomes the well known Wilson loop of size r x r (see 
e.g. Ref. [3T]). One can also fix the Coulomb gauge and 
omit the spatial Wilson line U(x, y;r). In this case one 
gets correlation function of two Wilson lines of extent r. 
The static meson correlation function has the following 
spectral representation 31j 



The coefficients c n depend on the choice of the operator, 
while the energy levels E n (r) corresponding to all possible 
states containing static quark and anti-quark do not [31 j . 
If truncated to include some number of terms, the above 
spectral representation can in principle be used to deter- 
mine the energy levels of QQ states in QCD at zero tem- 
perature. The lowest energy level of QQ pair is commonly 
referred to as the static potential, while the higher en- 
ergy levels are called the hybrid potentials. At sufficiently 
high temperatures, however, Eq. ([2]) is of little use since 
r < 1 /T and one has to consider infinite number of terms 
in Eq. As mentioned above at finite temperature the 
energy of a static QQ pair has also an imaginary part. It 
is not clear how this imaginary part would manifest itself 
in the Euclidean correlation function. It was suggested to 
generalize the spectral representation of the static meson 
correlators in terms of a temperature dependent spectral 
function [2"Tll2"g] 



G(r, 



dcoa r (oj, T)e 



(3) 



At zero temperature Eq. ([2]) and Eq. ([3]) are equivalent 
and the spectral function <j(uj, T) is just the sum of delta 
functions. The imaginary part of the static energy now 
is encoded in the width of the spectral function. While 
the spectral function ay(w,T) depends on the choice of 
the meson operator 0(x,y;r), the peak position and its 
width will not as long as there is a well defined bound 
state peak. Another way to arrive at Eq. is to consider 
the correlation function of static meson operators in real 
time 



D>(t) = Tr(O(t)O(0)exp{-l3H)) 

/— exp(— iujt)D > {to) 
-oo 2 7T 



(4) 



and then continue to imaginary time D > (t) — > D > {— ir). 
Using the standard definition of the spectral function 



D>(uj) - D<(m) 
2^r 



(5) 



G(r, t) = J~] Cn(r) exp(--E w (r)r), r 



\ x - V\ 



(2) 



and recalling that D < = in the case of static quarks (see 
discussions in Ref. [53]) we get Eq. ©. Note, that unlike 
quarkonium spectral function a r (uj) is not odd function 
of lu. From the point of view of MEM reconstruction of 
the spectral function oy(w,T) from Eq. ([3]) is not much 
different from the reconstruction of the quarkonium spec- 
tral functions. However, the structure of the spectral func- 
tion in the static case could be simpler and the cancella- 
tion of the temperature effects in the correlation function 
due to subtle interplay between the bound state peak and 
the continuum is not expected to happen. As a result we 
should see a more pronounced temperature dependence of 
the correlation functions. The above discussion holds for 
r < 1/T. The case t — 1/T should be considered sepa- 
rately because it is related to the free energy of static QQ 
pair. 

The correlation function evaluated at the maximal Eu- 
clidean time extent gives the so-called singlet free en- 
ergy F x {r,T) = —T\nG(r, r = 1/T). The singlet free 
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Fig. 1. Static meson correlators as function of r for two spatial separations: rT =1/2 (upper panels) and rT — 1 (lower panels). 
The numerical results are shown for two temperatures T = 147 MeV (left) and T = 266 MeV (right). 



energy has been studied in the past in pure gauge the- 
ory [^JIMIIMIISTIISSIIM] and in QCD [371IM1I3S1I1S] • Most 
of these studies use Coulomb gauge correlators but in 
Refs. [3"2ll3"rj] also Wilson loops have been considered. Un- 
like the true free energy of a static quark anti-quark pair 
F(r,T), defined in terms of Polyakov loop correlators, 
Fi (r, T) depends on the choice of the gauge or the choice 
of U(x,y;r) and thus is not physical. However, in the 
limit of the infinite separation it gives the QQ free en- 
(T) = linv-j.oo F(r, T) = lim^ Fx (r, T) . The 
singlet free energy is gauge independent in hard thermal 
loop (HTL) approximation [41] and coincides with the real 
part of the static energy in this approximation [241125] . 
The singlet free energy also arises in the context of short 
distance, rT <C 1 behavior of the Polyalov loop correlator 
and its decomposition into singlet and octet contribution 
within the pNRQCD framework [42]. Here it is defined 
as the singlet field correlator and thus is independent of 
U (x, y; t) or gauge fixing. It is equal to the zero temper- 
ature static potential up to power corrections, which are 
well defined and calculable at any order of perturbation 
theory. The singlet free energy defined in Coulomb gauge 
or in terms of rectangular (unsmeared) Wilson loops was 
studied using resummed HTL perturbation theory [43] . 
This calculation revealed that new type of temperature 
dependent divergences, so-called intersection divergences 



appear in the case of Wilson loops. These divergences 
have been shown to arise due to mixing of cyclic Wil- 
son loops with the Polyakov loop correlators [44] and can 
be removed if a proper renormalization procedure is im- 
plemented. This should make a quantitative comparison 
of the perturbative results with lattice data on the singlet 
free energy possible, thus providing a stringent test of the 
weakly interacting nature of QGP. 



3 Numerical results 

3.1 Lattice approach and parameters 

We calculated static meson correlation functions on 24 3 x 6 
lattice. Because of relatively small N T small statistical er- 
rors can be achieved. The gauge configurations used in 
our study have been generated by the HotQCD collabora- 
tion using a combination of the tree-level improved gauge 
action and Highly Improved Staggered Quark (HISQ) ac- 
tion [45 . This combination of the gauge action and quark 
action was referred to as the HISQ/tree action in Ref. [5S] 
but here we refer to it as the HISQ action for simplicity. 
The Highly Improved Staggered Quark action was first 
discussed in Ref. [13] . The gauge configurations have been 
generated using rational hybrid Monte-Carlo algorithm 



4 



A. Bazavov, P. Petreczky: Static meson correlators in 2+1 flavor QCD at non-zero temperature 



Table 1. Parameters of the lattice calculations 
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Fig. 2. The effective potential for rT — 1/2 as function of r 
at the lowest temperature T = 147 MeV. The open and filled 
diamonds correspond to Coulomb gauge correlators at zero and 
non-zero temperature, respectively. 



[47] . The algorithmic details of dynamical HISQ simula- 
tions can be found in Ref. 05] • The simulations have been 
performed for the physical value of the strange quark mass 
m s and light quark masses mi — m s /20. This light quark 
mass corresponds to the pion mass of 160 MeV in the 
continuum limit [45] . The parameters of the lattice simu- 
lations including the lattice gauge coupling f3 = 10/ g 2 and 
the strange quark mass in lattice units are shown in Table 
[1] along with the corresponding temperatures. As one can 
see from the table we consider a wide temperature range 
across the chiral crossover, which for N T — 6 occurs at 
T c = 171(1) MeV 05J. The last column of the table shows 
the accumulated statistics for each (3 value in terms of 
molecular dynamics time units (TU). Static meson corre- 
lators have been calculated after each 10 TUs. The lattice 
spacing a is determined from the r\ parameter defined in 
terms of the zero-temperature static potential as 



and we use the value ri = 0.3106 fm 05] , We use the pa- 
rameterization of the lattice spacing and the quark masses 
as functions of the gauge coupling (3 along the lines of con- 
stant physics that are given in Ref. [45] . 

Typically in lattice calculations a straight line path is 
used in the definition of the static meson operator in Eq. 
|T]), i.e. rectangular Wilson loops being calculated on the 
lattice. The main problem with calculating rectangular 
Wilson loops on the lattice and extracting potential from 
them is the large noise. To reduce the noise and improve 
the signal for the ground state potential smeared gauge 
links are used. The simplest example of smeared gauge 
links is the well known APE smearing [50]. In this case a 
3-link staple is added to each elementary lattice gauge link 
with some coefficient and the resulting sum is projected to 
SU(3) gauge group. Another smearing called hyper-cubic 
or HYP smearing was introduced in Ref. [5T] . In this case 
the smeared links are constructed from APE smeared links 
in lower dimensional subspace, such that the entire con- 
struction of the smeared links stays within the hypercube 



[51] . This smearing procedure depends on 3 parameters in 
4 dimensions, one for each level of APE smearing. Both 
APE and HYP smearing can be applied iteratively to in- 
crease the signal to noise ratio. Typically several steps of 
APE and HYP smearing are used. The HYP smearing 
turned out to be more efficient than the APE smearing in 
the calculations of the static potential in pure gauge the- 
ory 51. The APE smearing was used in calculating the 
static potential in 2+1 flavor QCD (see e.g. [52"ll53"] h while 
HYP smearing was used for example in Refs. [53]. An al- 
ternative approach to reduce the noise in the calculations 
of the static potential is to fix the Coulomb gauge (see 
e.g. 05|[551). Both approaches turned out to be equally 
efficient. 

In the present study we calculated static meson cor- 
relators G(t, r) at non-zero temperature using smeared 
spatial gauge links as well as Coulomb gauge. The smear- 
ing of the gauge links was restricted to spatial directions 
only, i.e. no temporal links enter the construction of the 
smeared links. We used several iterations of APE and HYP 
smearing. The coefficient of the staple for the APE smear- 
ing was c = 0.4, while for the HYP smearing we used the 
same smearing parameters as in Ref. [51], i.e. a-i = 0.3 
and «3 = 0.6 in the notation of Ref. [5T]. The parameter 
cti does not enter because we consider spatial directions 
only. We used 5, 10 and 20 iterations of APE smearing 
and 1, 2 and 5 iterations of HYP smearing. 

3.2 Time dependence of static correlatots 

As the first step toward understanding the temperature 
dependence of the static energy one should examine the r- 
dependence of static meson correlators for various temper- 
atures and different spatial separations r. As the medium 
effects on static meson correlators are expected to depend 
on how r and r compare to the inverse temperature we 
will scale these variables by the temperature when pre- 
senting our numerical results. Our numerical results for 
the static meson correlators as function of r with different 
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Fig. 3. The effective potential corresponding to various correlators at different temperatures for rT = 1/2 as function of r. The 
open and filled diamonds correspond to Coulomb gauge correlators at zero and non-zero temperature, respectively. 



smeared links as well as for Coulomb gauge are shown in 
Fig. [T]for two spatial separation rT — 1/2 and rT = 1. 
We also compared our results with zero temperature re- 
sults obtained in Coulomb gauge [35J. Note that unlike 
quarkonium correlators static meson correlators are not 
subject to periodic boundary conditions and thus can be 
studied also for tT > 1/2. As expected we see significant 
in-medium modification of the static meson correlators in 
contrast to quarkonium correlators. The modifications be- 
come larger with increasing r and increasing temperature: 
they are negligible at rT = 1/2, T = 147 MeV and are 
the largest for rT = 1 and T = 266 MeV. While the 
correlators in Coulomb gauge and smeared Wilson loops 
appear to be similar, unsmeared Wilson loops show quite 
different behavior. They appear to have more curvature 
as function of r and also have a different slope. One of the 
most prominent features of the finite temperature data is 
that the static meson correlator does not seem to follow 
an exponential decay at large r but flattens off or slightly 
increase around tT ~ 1. 

The increase around tT ~ 1 can be interpreted as a 
small contribution to the spectral function <j(u)) at neg- 
ative frequencies which arises due to periodic boundary 
conditions on the gluon background even-though static 
quarks propagate only forward in time, i.e. correspond to 
positive frequency. At least such interpretation seems to 



be supported by MEM analysis of Wilson loops [IS] and 
bottomonium correlators obtained in NRQCD [Tj5]. An- 
other way to understand the increase of the spatial meson 
correlator at tT = 1 is to inspect the HTL resummed per- 
turbative result for rectangular Wilson loop [23]. The re- 
summed 1-loop correction to the rectangular Wilson loop 
consist of three terms: a term linear in r that is propor- 
tional to the real part of the potential, a constant term and 
a term periodic in r [24] . After a proper re-exponentiation 
the term linear in r will correspond to the exponentially 
decaying part, while the two other terms can explain the 
behavior at tT ~ 1. 

While different correlators have different behavior in t 
they all converge at a similar value for tT = 1. For very 
large distances in fact this is expected, as — Tin G(r = 
1 /T, t) gives the free energy of static quark anti-quark 
pair and therefore all correlators should be equal (see the 
discussion in the previous section). The convergence of 
different correlators for rT — 1 is due to strong screening 
effects, i.e. the fact that — T In G(t = 1/T, r) is not very 
different from i 7 ^. 

At zero temperature the large r behavior of G(r, t) is 
dominated by the static potential. Therefore it is custom- 
ary to study the effective potential defined in terms of the 
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Fig. 4. The effective potential corresponding to various correlators at different temperatures for rT — 1 as function of r. The 
open and filled diamonds correspond to Coulomb gauge correlators at zero and non-zero temperature, respectively. 



following ratio 

aV eff (r, t) 



In 



G(r,r/q) 
G(r,r/a + iy 



(7) 



Here a denotes the lattice spacing. The effective potential 
is expected to decrease with increasing r, eventually reach- 
ing a plateau that corresponds to the value of the static 
potential. How quickly V e f / reaches a plateau depends on 
the choice of the static meson operator in Eq. ([1]), using 
more smearing in the gauge connection U(x^y;r) typi- 
cally suppresses the contribution from the excited states. 
Using the Coulomb gauge is also an effective method for 
suppressing the excited state contributions. The effective 
potential turns out to be a useful tool to examine the 
behavior of the static meson correlators in more detail. 
In Fig. [2] we show our numerical results for V e f t calcu- 
lated for different smearing of the spatial links as well 
as for Coulomb gauge for the lowest temperature of 147 
MeV and distance rT = 1/2. We do not show V e ff from 
unsmeared Wilson loops as the corresponding numerical 
results are too noisy. We compare our numerical results 
with the zero temperature results in Coulomb gauge [4"5] 
which approach a plateau within the considered r range. 
At the lowest temperature we do not expect large mod- 
ification of the static potential. Indeed the temperature 
effects on V e ff are small for tT < 1/2 and indicate a 



downward shift of the static potential of about 10 MeV. 
This is consistent with the study of the static potential 
in pion matter, where a similar downward shift has been 
observed [SI] (see also Ref. [S3] for a related study for 
quarkonium). The figure also shows that Coulomb gauge 
results and smeared Wilson loops give similar results if 
sufficient number of smearings steps is used, i.e. 5 or more 
smearing steps in the case of APE smearings and more 
than two steps for HYP smearings. For one or two steps 
of HYP smearings the excited state contribution, seen as 
the curvature at small r, is significant. While the effective 
potential reaches a plateau at zero temperature within the 
studied r-range, this is not the case for the lowest tem- 
perature. The decrease in V e ff at rT = 2/3 is particularly 
striking and difficult to explain in terms of the expected 
in-medium modifications of the static energy. Similar re- 
sults have been obtained for T = 155 MeV, 162 MeV and 
178 MeV. 

Now let us discuss the behavior of V e f / above the chi- 
ral transition temperature. We show our numerical results 
for rT = 1/2 and rT = 1 in Fig. [3] and Fig. H respec- 
tively. The features of V e // discussed above also hold at 
higher temperatures, in particular, V e ff does not show a 
plateau and there is a rapid decrease around tT — 2/3. 
The effective potential above T c is always smaller than the 
T = result and the deviations increase with increasing 
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Fig. 5. The singlet free energy as function of distance r for 
different temperatures T < 200 MeV calculated in Coulomb 
gauge (top) and using smeared Wilson loops with 5 steps of 
HYP smearing (bottom). We also show the zero temperature 
static potential for /3 = 6.664 45 . 



a simple form CFa s m,D, while it vanishes for rT <C 1. For 
the relevant temperature range a s nir) ~ 1. Therefore in 
our modeling efforts we assume that IvaV — T for rT = 1 
and ImV = T/2 for rT = 1/2. The value of IW was 
adjusted to match the value of V e // at rT = 1/2. The 
corresponding results are shown in Figs. [3J and 2] as solid 
black lines. For rT = 1 we also allow for a larger value of 
the imaginary part of the potential, namely ImV = 3T. 
The corresponding results are shown as the dashed lines 
in Figs. [3] and [4] As one can see from the figures the r 
dependence of V e ff in general is not described well by 
the above form of the spectral function. For rT = 1/2 
the Breit-Wigner form captures some features of the r 
dependence but for rT = 1 the observed curvature of V e / / 
is much larger than the one that can be provided by the 
Breit-Wigner form, even if we assume that ImV = 3T. The 
Breit-Wigner form describes the spectral function well in 
the vicinity of the pole. For large values of InV the Breit- 
Wigner form has long tails that extend to the u regions 
where the form is not valid, e.g. Eq. ([5]) obtained from 
pNQRCD is not valid for large values of ui. Therefore we 
also consider a Gaussian Ansatz with the same width as 
used in the Breit-Wigner form and calculate G(r,r). In 
this case we get very small T-dependence of V e / / . Thus the 
r dependence seen in Figs. [3J and 2] for the Breit-Wigner 
form is generated by its tails. Therefore it seems unlikely 
that the T-dependence of V e / / is due to the imaginary part 
of the static potential, instead there could be an additional 
contribution to the static correlator that arises from the 
periodic boundary. 

As discussed above, based on HTL resummed results 
we expect that Wilson loops at high temperatures will con- 
sist of two contributions: one exponentially decaying in r 
and another one with more complex r behavior. Therefore 
we tried to fit the r dependence of the static correlators 
by a simple form 



temperatures and distance suggesting that the static po- 
tential decreases with increasing temperature. Note that 
we do not present our results for V// from unsmeared 
Wilson loops as these show much stronger T-dependence 
and would change the scale of Figs. [3J and 0] substantially. 

At temperatures significantly higher than the transi- 
tion temperature we expect the imaginary part of the po- 
tential to play an important role. Therefore it is natural 
to ask whether the above features of V e // are due to the 
imaginary part. According to thermal pNRQCD the spec- 



Aexp(-Vr) + cexp(-A(l/T - t)) 



(9) 



tral function defined in terms of D' 



i.e. a r 



has the Breit-Wigner form [33] 

<T r (uj) = - 



D>{u)/{2k) 



ImV 



7r (w - ReV) 2 + ImV 2 



(8) 



(c.f. Eq. (96) in Ref. [23]). To model the effects of ImV on 
the effective potential we calculated the static correlators 
using the above form for the spectral function and Eq. 
(J3J together with some plausible values of ImV inspired 
by weak coupling calculations. The r-dependence of the 
imaginary part of the potential is quite complicated even 
in the weak coupling regime. For large r, however, it has 



for 1/2 < t < 2/3. The fit form works well but the ex- 
tracted value of V depends somewhat on the chosen fit 
range. However, it is always smaller than the T — static 
potential. This together withthe analysis of V e ff suggests 
that the energy of static QQ pair is decreasing with in- 
creasing temperature contrary to the results obtained in 
Ref. [28] based on calculations in SU(3) gauge theory and 
MEM. We also find that the non-exponential contribution 
is relatively small. Namely we find that c is about 1000 
times smaller than A, however, it increases with increasing 
temperature. 



3.3 Singlet free energy 

As discussed in section 2 F\{r,T) = —TlnG(r,T) gives 
the singlet free energy of static QQ. We have calculated 
the singlet free energy F\ using Coulomb gauge and smeared 
as well as unsmeared Wilson loops. We start the discussion 
of our numerical results first considering the low temper- 
ature region T < 200 MeV. This region of course also 
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contains temperature values above the chiral crossover. 
Our results for Coulomb gauge and smeared Wilson loops 
are shown in Fig.[5Jand compared to the zero temperature 
static energy calculated in Ref. [3H] at (3 — 6.664. We do 
not show the results from unsmeared Wilson loops because 
they are too noisy. It is suffice to say that they are above 
the Coulomb gauge result. All the smeared Wilson loops 
give very similar results for F±, i.e. the details of smearing 
do not seem to matter. The singlet free energy does not 
show a significant change across the chiral crossover and 
the temperature effects are relatively small for r < 0.65 
fm. In the case of Coulomb gauge Fi is slightly larger than 
the zero temperature static potential at intermediate dis- 
tances 0.3 fm < r < 0.8 fm. This is qualitatively similar to 
the findings in pure gauge theory [54"ll551l5"Tj . As discussed 
in Ref. [55J this is due to the fact that c\ < 1. In the case 
of smeared Wilson loops F\ is always smaller than the zero 
temperature static energy since c\ is expected to be close 
to unity [55] . 

Our results for the high temperature region, T > 200 
MeV are shown in Fig. [5J for Coulomb gauge and un- 
smeared Wilson loops. As expected at very small dis- 
tances F\ is temperature independent and coincides with 
the T = potential. For T > 200 MeV the numerical 
results for unsmeared Wilson loops have quite small sta- 
tistical errors. Therefore one can use our numerical results 
to test the perturbative predictions for cyclic Wilson loops 
[431144] . Note that for Coulomb gauge the singlet free en- 
ergy does not overshoot the zero temperature static poten- 
tial similarly to the case of pure gauge theory [5411571156] . 
This may imply that the coefficients c„ are not important 
for the temperature dependence of F\ at high tempera- 
tures. As the temperature increases the singlet free energy 
reaches the plateau at smaller and smaller distances con- 
sistent with expectations based on color screening. Figure 
[6] shows that F\ obtained from unsmeared Wilson loops is 
systematically larger than F\ obtained in Coulomb gauge 
at intermediate distances. This could be due to the inter- 
section divergences of cyclic Wilson loops [44]. Smeared 
Wilson loops on the other hand give results which are 
slightly smaller than Coulomb gauge results at interme- 
diate distances, except for one iteration of HYP smear- 
ing. There is no dependence on the smearing if more than 
one smearing steps are used. To make this point clear in 
Fig. [7] we compare F\ obtained for different correlators at 
T = 281 MeV. Similar results have been obtained at other 
temperatures. 

To explore the large distance behavior of the singlet 
free energy we consider the following combination 



5 1 (r,T) = -r(F 1 (r,T)-F 00 (T)). 



(10) 



At leading order this combination should decay exponen- 
tially ans its exponential decay is governed by the lead- 
ing order Debye mass = gT-^3/2. The perturbative 
analysis of F\ quantity at next-to-leading order was first 
performed in Ref. [35] and in the case of Colomb gauge 
revealed a power-law behavior. However, as mentioned al- 
ready the renormalization of the static meson correlators 
at t = 1/T is non-trivial. We find that Si(r,T) decays 
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exponentially for r > 1/T. The corresponding screening 
mass m e /T shows only mild temperature dependence for 
T > 300 MeV consistent with expectations. Furthermore, 
we find that ttie/T for unsmeared Wilson loops is sig- 
nificantly larger than for Coulomb gauge. The screening 
masses obtained from smeared Wilson loops are similar or 
slightly smaller than the ones obtained in Coulomb gauge. 
The only exception is the case when one iteration of HYP 
smearing is used, where to^/T is closer to the the value 
obtained from unsmeared Wilson loops. 



4 Conclusions 

In this paper we studied various static meson correlators 
at non-zero temperature. We analyzed the temperature 
dependence of these correlators as well as their depen- 
dence on the Euclidean time r. The r-dependence of the 
correlators has been analyzed in terms of the effective po- 
tential V e ff for various spatial separations r. For large 
Euclidean times V e t / gives the static potential in the zero 
temperature case. We found that static meson correlators 
defined in Coulomb gauge and in terms of smeared Wilson 
loops show quite similar behavior. Contrary to the zero 
temperature case the effective potential does not reach a 
plateau in the studied t window even at the lowest tem- 
peratures. We explored different possibilities to explain 
this behavior including the presence of an imaginary part 
of the static energy. We conclude that this behavior is due 
to the additional contribution not related to the static po- 
tential that arises from the periodic boundary condition. 
To extract the potential at non-zero temperature we will 
need to extend this calculation to larger N T . We also cal- 
culated the singlet free energy and found that its behavior 
is qualitatively similar to the previous findings obtained in 
pure gauge theory. Comparing our numerical results with 
perturbative prediction on F\ 42,43,44 will be useful to 
clarify the strongly or weakly interacting nature of the 
QGP. 
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